Andreev-Klein reflection in graphene ferromagnet-superconductor junctions 
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We show that Andreev reflection in a junction between ferromagnetic (F) and superconducting 
(S) graphene regions is fundamentaUy different from the common FS junctions. For a weakly doped 
F graphene with an exchange field h larger than its Fermi energy Ep , Andreev reflection of massless 
Dirac fermions is associated with a Klein tunneling through an exchange field p-n barrier between 
two spin-split conduction and valence subbands. We find that this Andreev-Klein process results 
in an enhancement of the subgap conductance of a graphene FS junction by h up to the point at 
which the conductance at low voltages eV ^ A is greater than its value for the corresponding non- 
ferromagnetic junction. We also demonstrate that the Andreev reflection can be of retro or specular 
types in both convergent and divergent ways with the reflection direction aligned, respectively, closer 
to and farther from the normal to the junction as compared to the incidence direction. 

PACS numbers: 74.45.-l-c, 73.23.-b, 85.75.-d, 74.78. Na 
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Transmission of low energy electrons through a normal- 
metal -superconductor (NS) junction is realized via a pe- 
culiar scattering process, known as Andreev reflection 
(AR) 1]. In AR an electron excitation with energy e 
and spin direction a upon hitting the NS-interface is con- 
verted into a hole excitation with the same energy but 
opposite spin direction a — —a. Under AR the momen- 
tum change is of order e/v-p which is negligibly small for 
a degenerate N metal with large Fermi energy E-p ^ A. 
Thus the hole velocity is almost opposite to the veloc- 
ity of the incident electron (since a hole moves opposite 
to its momentum), which implies that Andreev process 
is retro reflection. Andreev reflection results in a finite 
conductance of a NS junction at the voltages below the 
superconducting gap A [3|. 

The fact that the Andreev reflected electron-hole be- 
long to different spin-subbands has an important conse- 
quence for Andreev conductance when N metal is a ferro- 
magnet (F) . The exchange splitting energy h of F-metal 
induces an extra momentum change 2/i/uf of the re- 
flected hole which diminishes the amplitude of AR. As the 
result the subgap Andreev condutance of ferromagnet- 
superconductor (FS) junctions decreases with increasing 
h and vanishes for a half-metal F with h = Ep [3[ . Sup- 
pression of AR at FS interface is a manifestation of the 
common fact that ferromagnetism and spin singlet su- 
perconductivity are opposing phenomena. In this letter, 
however, we show that the situation differs significantly 
if the FS junction is realized in graphene, the recently 
discovered two-dimensional (2D) carbon atoms arranged 
in hexagonal lattice 4, 5,, 6]. We find that in a graphene 
FS junction, the exchange interaction can enhance the 
subgap Andreev conductance, depending on the doping 
of F graphene. In particular we show that at low voltage 
eV <§:_ A the conductance of a graphene FS junction with 
a strong exchange field h ^ Ep is larger than its value 
for the corresponding NS structure. We explain this ef- 
fect in terms of Andreev-Klein reflection in which the 
superconducting electron-hole conversion at FS interface 



is accompanied with a pseudo-relativistic Klein transmis- 
sion through an exchange built p-n barrier between the 
two spin-split conduction and valence subbands. 

Graphene is a zero-gap semiconductor with its con- 
ical valnce and conduction bands touching each other 
at the corners of hexagonal first Brillouin zone, known 
as Dirac points. The carrier type, (electron-like (n) or 
hole-like (p)) and its density can be tuned by means of 
electrical gate or doping of underlying substrate. Due 
to the connection between this specific band structure 
and the pseudo-spin aspect which characterizes the rela- 
tive amplitude of electron wave function in two trigonal 
sublattices of the hexagonal structure, the charge carri- 
ers in graphene behave like 2D massless Dirac fermions 
with a pseudo-relativistic chiral property [a, Ifl lD, Is], Isl] . 
Currently intriguing properties of graphene, which arise 
from such a Dirac like spectrum, have been the subject 
of intense studies 10|, [ill, |l2 1 . 



Among others, peculiarity of Andreev reflection in 
graphene NS junctions was predicted by Beenakker 
13l . Il4| . S region with high carrier density can be pro- 



duced by depositing superconducting electrodes on top 
of a graphene sheet T5j|. It was demonstrated that un- 
like the highly doped graphene or a degenerate N metal, 
for undoped graphene with A ^ Ep the change in the 
momentum upon AR could be of order of the momentum 
of the incident electron. In this limit the dominant pro- 
cess is AR of an electron from the conduction band into 
a hole in the valence band in which the reflection angle 
(versus the normal to the NS interface) is inverted with 
respect to the incidence angle, making Andreev process 
a specular reflection. Transition from retro reflection at 
A ^ i?F to specular AR at A ^ ii^p is associated with 
an inversion of the voltage dependence of the subgap An- 
dreev conductance. 

Recently proximity induced ferromagnetism was exper- 
imentally realized in graphene spin- valve structures [l6| . 
Intrinsic ferromagnetic correlations were also predicted to 
exist in graphene sheets 17| and nanoribbons 1§\. For a 



pure F graphene sheet the exchange energy shifts the nor- 
mal Fermi level at Dirac point {E-p = 0) upward (down- 
ward) by h in its spin-up (down) subband. A important 
consequence of the gapless Dirac spectrum is that this 
shift makes the up and down spin carriers to be electron- 
like (n type) and hole-like (p type), respectively [17 1. 
Concerning the transport between up and down spin- 
subbands the exchange field, thus, operates as a p-n po- 
tential barrier. The similar situation happens for a doped 
F graphene {Ep ^ 0) samples with large exchange ener- 
gies h :s> i?F(see the right inset in Fig. [2). At a graphene 
FS interface AR converting electron-hole from different 
spin subbands will bring this exchange correlations built 
p-n barrier into effect. Already reflectionless transmis- 
sion of chiral electrons through wide and high normal 
graphene p-n barriers was demonstrated 19J, |20|, |2l|, \2^ . 
This effect called Klein tunneling is analogous to the cor- 
responding effect in quantum relativistic theory [23| . We 
show that the spin Klein tunneling at graphene FS junc- 
tion leads to an enhancement of the amplitude of AR and 
the resulting Andreev conductance by the exchange field. 
This finding specific to graphene is in striking contrast to 
the behaviour of Andreev conductance of a FS junction 
in the ordinary metals. 

Concerning the connection between the incidence and 
reflection directions, we further demonstrate variety of 
Andreev processes taking place in graphene FS junctions. 
For an incident spin-a electron with the velocity direc- 
tion angle 0^ versus the normal to the junction, Andreev 
reflection can be of retro or specular types indicated, 
respectively, with or without a sign change in the re- 
flection angle, sign[(/()^] = ±sign[(/>o-], in both convergent 
10^1 < |(/)cr| and divergent |(/)^| > |(/)cr| ways. The type of 
AR depends on the ratio h/ E-p, the spin a and energy e 
of the incident electron. 

We consider a wide graphene FS junction normal to 
X-axis with ferromagnetic region for x < and highly 
doped superconducting region for x > 0. In the F region 
the two up and down (cr = ±) spin subbands are split by 
2h, such that the spin a excitation spectrum versus 2D 
wave vector k^ = {ka-,qa-) is give by 



Ea = \Ef ± hvlkal + ah\, 



(1) 



where the two branches ± of the spectrum originate 
from the valence and conduction bands, respectively. 
In S region h — and the superconducting correla- 
tions are characterized by the order parameter A which 
is taken to be real and constant. For calculation we 
adopt Dirac-Bogoliubov-de Gennes equation [l3| which 
describes the superconducting correlation between mass- 
less Dirac fermions with different valley indices. In the 
presence of an exchange interaction it has the form 
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where Hq = —ihvY{axdx+Uydy) — U{r)—EF is the Dirac 
Hamiltonian and the potential energy U{r) = Uq ^ Ep 
in S and U{r) — in F; <Tx and Uy are Pauli matrices in 
the pseudospin space of the sublattices. 

Within the scattering formalism we find the spin- 
dependent amplitude of Andreev and normal reflections 
from FS interface. An incident spin-u electron from 
left to FS interface with a subgap energy e < A and 
the incidence angle (pa- — a,rcsm{hvFq/{e + Ep + ah)) 
can be either normally reflected or Andreev reflected 
as a hole with opposite spin a along the reflection an- 
gle (f)'g. = a,rcsin{hvFq/{e — Ep + ah)). Denoting the 
amplitude of normal and Andreev reflections, respec- 
tively, by To- and rAa the wave function inside F is writ- 
ten as *F(T == V'+o- + ^crfc + TAcriJ-'ha, where ipf^ « 
e:icp{iq±ikax) x (exp(=FJ(/)£r),±l,0, 0) and ipj^- cc exp(zg — 
ikgx) X (O,O,exp(i0^), 1), are the eigenstates of Hamil- 
tonian ([2]) in F. The transmitted part of the electron 
into S, ^5cr = datps + baipg consists of two super- 
conducting quasiparticles whose wave functions ips ~ 
e'iq){iq+iks±x) x (exp(±i/?), ± exp(±i/3), 1, ±1) decay ex- 
ponentially as a function of x (/3 = arccos(e/A)). Match- 
ing the wave functions in F and S at the interface x = 
we obtain 



sec 
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{<P'-„ - <P„)I2 -f ^ tan/3cos(<?!)^ + 0^)/2 ' 
m{(p's + (pa)l2 -t- ztan/3sin(0;, - 0^)/2 
i{^'a-<Pa)l1 + i\.^.'i^(icQs{<P'a + <Pa)l2 ■ 



(3) 
(4) 



From the conservation of the j/-component wave vector 
q under the scattering we obtain the following relation 
between the incidence and reflection angles 
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e^Ep + ah 
e — Ef + ah' 



(5) 



(2) 



Inspection of the results given by Eqs. ([T|) and ([3])- 
([5]) reveals variety of spin-dependent Andreev processes 
in the FS graphene junction. Consider an spin-up elec- 
tron with the energy e above the Fermi level in the 
conduction band {Ep > 0) hitting the junction in an 
angle 0+ versus the normal to the junction (negative 
X-axis). It can be Andreev reflected as a hole with 
the same energy below the Fermi level in down-spin 
subband. From relation ([T]) we see that as long as 
e + h < E-p the hole is in the conduction band too, and 
AR is retro. Upon AR q and e are conserved but the 
magnitude of the momentum is changed by 2{e -\- K)/vf 
(see Eq. ([1])). These conditions result in the relation 

k'_ = J{e + h- £'f)2 - {e + h + Ep)'^ sin^ (/>+, for the x- 
component of the down spin hole which shows that above 
the critical angle 0!^ = arcsin(|e — Ep + h\/\e + E-p -^ h\), 
the hole wave function is evanescent and amplitude of 
AR vanishes. When e + h ^ 0, (p'^ = 7r/2 and the x- 
component of the reflected hole is the same as the inci- 
dent electron, so 6'_ = —6+. This deflnes a line in the 
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FIG. 1: Map of the Andreev reflection angle (?!)'_ for an up- 
spin electron incident with the angle 4>+ on the grapgene FS 
junction. It shows dependence on the electron energy e/ E-p 
and the exchange energy h/Ep in the scale of Fermi energy. 
In the shaded region between the two lines e/E-p + h/Ep — 
±1 the Andreev reflection is retro while in the white region 
it is specular. The line e/Ep + h/Ep — 0, on which the 
reflection is perfectly retro, divides the regions of convergent 
reflection(below) and divergent reflection(above). 



phase space of e/Ep and h/Ep (see Fig. [T]) on which 
AR is perfectly (incidence and reflection directions are 
aligned) retro. For a fixed e, increasing h results in an 
increase in 1 0'_ | , implying a divergent retro reflection (see 
Eq. ([5])). As h approaches Ep — e the angle (j>'_ -^ — 7r/2 
and simultaneously 0^ decreases monotonically to ap- 
proach zero. This implies that close to the line e+h ~ Ep 
only normally incident spin-up electrons have a finite AR 
amplitude. Thus in h/ Ep — e/ Ep phase diagram shown 
in Fig. [T] the region between two lines e + h = Q and 
e + h = Ep defines the phase of divergent retro AR. In- 
creasing h further we cross the line e + h = Ep which 
is associated with a sign change in the reflection angle 
from — 7r/2 to 7r/2 implying a transition to the regime of 
specular AR. For higher /i, <f>'_ decreases toward 0-|_ and 
(jfj^ increases. In the limit oi h + e '^ Ep, (j)'_ — > (j>^ and 
AR becomes perfectly specular. This defines the region 
of divergent specular AR above the line e + h = Ep . 

In a similar way we observe that starting from the per- 
fectly retro AR line and decreasing h for a fixed e, AR 
remains retro but becomes convergent with \4>'_\ < \4'+\- 
The retro convergent region —(f>+ < 4>'_ < Q extends be- 
tween the two lines of e + /i = and e + h ~ —Ep. Ap- 
proaching the line e + h — ~Ep, (/)'_—> 0. Upon crossing 
this line we again will have a retro to specular transi- 
tion, but this time in a continues way at </)'_ — 0. Below 
the line of <j>'_ = 0, is the region of convergent specular 
AR with < </)'_ <</)+■ For all convergent reflection 
region 0!j_ = tt/2. Note that at zero energy the retro or 
specular (depending on h/ Ep) AR will transform from 
divergent to convergent or vice versa upon changing the 



FIG. 2: Dependence of the Andreev conductance of a 
graphene FS contact on the exchange field h/Ep in units of 
the Fermi energy for a highly doped F, Ep 2> A, and at three 
bias voltages el//A = 0, l/\/2, 1 (solid, dashed and dotted 
curves, respectively). For h > Ep the conductance increases 
with h/Ep. Note that for an undoped graphene [Ep = 0) 
the conductance for each of the subgap voltages takes an ex- 
change field-independent value determined by its value for a 
doped sample in the limit h/Ep ^ 1. The insets show the 
configuration (being electron- like (n), hole- like (p) or neutral 
(D) at Dirac point) of the up and down spin-subbands for 
three different cases of h/Ep < 1, h/Ep = 1 and h/Ep > 1. 



spin direction of the incident electron (see Fig. [T|) . 

We have calculated the Andreev conductance of FS 
junction at zero temperature via Blonder-Tinkham- 
Klapwijk formula 
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d0<,cos0^(l-|r^p-|-|rAo-n, (6) 



where the spin-u normal-state conductance G^ ~ 
{2e^/h)N„{eV) and the density of states N„{e) = \e + 
Ep A- ah\W / {TThvp) iW is width of the junction). 

Dependence of the resulting Andreev conductance 
G/Gp {Gp = G+ + G-) on h/Ep is presented in Fig. 
[2]for three different bias voltages eV/A — 0, 1/a/2, 1 and 
for highly doped F graphene Ep/A 3> 1. For h/Ep < 1 
the conductance decreases monotonically below the cor- 
responding NS value G/Gp{h = 0) with h/Ep and van- 
ishes &t h — Ep. In this regime the up and down spin 
subbands are of the same n type (left inset of Fig. [2|) 
and the effect of exchange field is to impose a normal 
barrier against AR. The resulting momentum difference 
of Andreev reflected electron-hole diminishes AR ampli- 
tude and hence the Andreev conductance. M h — Ep 
the down spin subband is at Dirac point (middle inset of 
Fig. [2]) with a vanishing density of states which results 
in G/Gp = 0. For h > Ep Fermi level of spin down elec- 
trons is transferred into the valence band and thus the 
exchange field barrier finds a p-n characteristic with the 
height 2/i— _EF(right inset of Fig. [5]). In this case the con- 



ductance G/Gp increases monotonically with h/E-p- In 
the hmit of h/E-p ^ 1 the exchange barrier transforms 
to an almost perfect p-n barrier with the height ~ 2h 
which makes possible perfect transmission of chiral Dirac 
fermions between the two subbands. The enhancing An- 
dreev conductance reaches a limiting maximum which 
depends on the bias voltage. Importantly we see that 
at eF/A = this hmiting value G/Gp{h/Ep > 1) = 2 
is larger than the value for corresponding NS structure 
G/Gp{h/Ep < 1) = 4/3 as shown in Fig. H This ef- 
fect is in contrast to the common view that that ferro- 
magnetic ordering and the singlet superconductivity are 
exclusive phenomena. 

Finally, we note that with the recent successful induc- 
tion of superconductivity [l5| and spin-polarization [l6| 
in graphcnc sheets, realizing graphene FS junctions seems 
feasible. This will allow for the experimental observation 
of the above predicted effect. 

In conclusion, we have demonstrated the unusual fea- 
tures of Andreev reflection in graphene ferromagnet- 
superconductor junctions. We have shown that depend- 
ing on the ratio of the exchange field and Fermi energy 
h/Ep in F graphene and the energy and spin direction 
of the incident electron, the Andreev reflection may be 
of retro or specular types with possibility of being con- 
vergent as well as divergent. More fundamentally and 
in contrast to the common view, we have found that for 
h > Ep the exchange field enhances the Andreev conduc- 
tance of FS junction to reach a maximum value which at 
low bias voltages eV <C A is greater that its value for 
the corresponding non-ferromagnetic structure. We have 
explained this effect in terms of Andreev-Klien reflection 
in which the superconducting electron-hole conversion at 
FS interface is associated with a Klein tunneling through 
the exchange built p-n barrier between the spin-split con- 
duction and valence subbands. 



[9; 
[lo; 

[11 
[12; 

[is; 

[14 

[is; 
[16; 
[17; 

[is; 

[19 
[20 
[21 
[22; 
[23; 



A. F. Andreev, Sov. Phys. JETP 19, 1228 (1964). 
G. E. Blonder, M. Tinkham, and T. M. Klapwijk, Phys. 
Rev. B 25, 4515 (1982). 

M. J. M. de Jong and C. W. J. Beenakker, Phys. Rev. 
Lett. 74, 1657 (1995). 

K. S. Novoselov et ai, Sience 306, 666 (2004). 
K. S. Novoselov et ai, Nature 438, 197 (2005). 
Y. Zhang, Y. W. Tan, H. L. Stormer, and P. Kim, Nature 
438, 201 (2005). 

P. R. Wallace, Phys. Rev. 71, 622 (1947). 
J. C. Slonczewski and P. R. Weiss, Phys. Rev. 109, 272 
(1958). 

F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 
A. H. Castro Neto, F. Gui nea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim,'arXiv:0709.1163 (2007). 
A. K. Geim and K. S. Novoselov, Nature Mat. 6, 183 
(2007). 

M. L Katsnelson and K. S. Novoselov, Solid State Comm. 
143, 3 (2007). 

C. W. J. Beenakker, Phys. Rev. Lett. 97, 067007 (2006). 
C. W. J. Beenakker, arXiv:0710.3848 (2007). 
H. B. Heersche et ai, Nature 446, 56 (2007). 
N. Tombros et al, Nature 448, 571 (2007). 
N. M. R. Peres, M. A. N. Araiijo, and D. Bozi, Phys. 
Rev. B 70, 195122 (2004); N. M. R. Peres, F. Guinea, 
and A. H. Castro Neto, Phys. Rev. B 72, 174406 (2005). 
Y. -W. Son, M. L. Cohen, and S. G. Louie, Nature 444, 
347 (2006). 

T. Ando, T. Nakanishi, and R. Saito, J. Phys. Soc. Japan 
67, 2857 (1998). 

V. V. Cheianov and V. I. Fal'ko, Phys. Rev. B 74, 
041403(R) (2006). 

M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Na- 
ture Phys. 2, 620 (2006). 

C. W. J. Beenakker, A. R. Akhmerov, P. Recher, and J. 
Tworzydlo, Phys. Rev. B 77, 075409 (2008). 
O. Klein, Z. Phys. 53, 157 (1929). 



